Origin of chaos in soft interactions and signatures of non-ergodicity 
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The emergence of chaotic motion is discussed for hard, point-like, and soft colhsions between 
two particles in a one-dimensional box. It is known that ergodicity may be obtained in point-like 
collisions for specific mass ratios 7 = ni^jriix of the two particles, and that Lyapunov exponents are 
zero. However, if a Yukawa interaction between the particles is introduced, we show analytically 
that positive Lyapunov exponents are generated due to double collisions close to the walls. While 
the largest finite-time Lyapunov exponent changes smoothly with 7, the number of occurencies 
of the most probable one, extracted from the distribution of finite-time Lyapunov exponents over 
initial conditions, reveals details about the phase space dynamics. In particular the influence of the 
integrable and pseudointegrable dynamics without Yukawa interaction for specific mass ratios can 
be clearly identified and demonstrates the sensitivity of the finite-time Lyapunov exponents as a 
phase space probe. Being not restricted to two-dimensional problems such as Poincare sections, the 
number of occurencies of the most probable Lyapunov exponents suggest itself as a suitable tool 
to characterize phase space dynamics in higher dimensions. This is shown for the problem of two 
interacting particles in a circular billiard. 
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I. INTRODUCTION 

The investigation of the origin of chaotic motion in 
standard billiard models (like the Sinai billiard the 
Bunimovich stadium or the Annular billiard [31), has 
played a pioneering role since the very beginning of chaos 
theory. Usually in such models, ballistic chaotic mo- 
tion (single particle dynamics) is a consequence of the 
spatial billiard geometry. For interacting many-particles 
systems, which appear in many areas of physics, chaotic 
motion can be generated from the combined effect of ex- 
ternal forces and mutual interactions. In order to un- 
derstand how chaotic motion emerges as a consequence 
of the interaction between particles, a simple billiard, 
namely a one-dimensional box, will be used. Since in this 
case the boundary alone cannot induces irregular motion 
of the two particles inside, the role of the interaction for 
the generation of chaotic motion becomes clear. 

Interacting particles inside billiards can be used to 
model electrons in quantum dots. Electrons are con- 
fined inside a disk and are affected by the surrounding 
material which composes the semiconductor Q. In fact, 
the composition of the surrounding material may destroy 
the long Coulomb repulsion between electrons and also 
change the effective mass between particles The in- 
fluence of both effects on quantum energy levels and/or 
electrons dynamics it not obvious. However, they can be 
studied in detail for a physical model where kinetic and 
potential energy of the particles can be varied indepen- 
dently. This is achieved by a parameter which controls 
the range of interaction between the particles (electrons) 
and by varying the mass ratio 7 = 7712/7711 of the parti- 
cles. 

In this paper we will study the classical dynamics of 



two interacting particles inside a one-dimensional billiard 
as a function of 7 and as a function of interaction range 
between particles. A Yukawa interaction between parti- 
cles is assumed. Such a system has been considered clas- 
sically and quantum mechanically Q for the case of 
equal masses. In order to calculate the spectrum of Lya- 
punov Exponents (LEs), the dynamics in tangent space is 
determined explicitly. In the limit of very short range of 
interaction, this system should approach the hard-point 
collision case analyzed originally by Casati et al Q . De- 
spite ergodic dynamics for the case of point-like colli- 
sion at specific 7, vanishing LEs are a consequence of the 
linear instability of this system Such linear unsta- 
ble systems have become a topical problem in statistical 
mechanics (see also about the origin of diffusion in 
non-chaotic systems [HI). 

In the case of a Yukawa interaction, the repulsion be- 
tween particles at collisions with the walls are shown to 
generate positive LEs. This is shown explicitly by deter- 
mining the dynamics in tangent space. While the mean 
value of the largest finite-time LE only quantifies the 
degree of chaoticity of the system, the number of oc- 
curencies of the most probable LE, extracted from the 
distribution over initial conditions, is shown to give sig- 
nificant information about regular structures and sticky 
(or trapped) [l2| trajectories in phase space. This distri- 
bution is determined numerically as a function of 7. 

The plan of the paper is as follows. Section [Tl] re- 
views the main results from the problem of two hard 
point particles in a one-dimensional box. In section IIIII 
chaotic motion emerges with the introduction of the soft 
Yukawa interaction between the particles. Positive LEs 
can be generated from analytical expressions obtained for 
the dynamics in tangent space. The distribution of the 
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largest finite-time LE calculated over the phase space of 
initial conditions is used to reveal the underlying dynam- 
ics. Section |lll] discusses the distribution of the largest 
finite-time LE for the case of two interacting particles in 
a circular billiard. The paper ends with conclusions in 
section |Vl 



II. TWO PARTICLES IN A ID-BOX WITH 
HARD POINT-LIKE COLLISIONS 

Two particles in a ID-box with hard point-like colli- 
sions, also called two-particle hard point gas, can also 
be treated as a particular case of the motion of three 
particles on a finite ring which can be mapped 

onto the motion of a particle in a triangle billiard [15| . 
In such systems the Lyapunov exponent is zero Q and 
the whole dynamics can be monitored by changing the 
angles of the triangle billiard. These angles are functions 
of the masses ratio between particles. Such triangle bil- 
liards have also been applied to study energy diffusion in 
one-dimensional systems 0] . Although the connection to 
the triangle billiard is very usefuU to gain insights about 
collision properties of the particles, it is not needed for 
the purpose of the present work. 

As follows we summarize the main results obtained by 
Casati Q, which are similar to those observed in the 
triangle billiard. Using Poincarc sections they 3 showed 
that the dynamics is non-crgodic if is a rational multiple 
of TT, where 



genus g ~ 2, the possible values of 9 are \li 



•TT, ^TT 



and 



the mass ratios are 7 ~ 0.1 and 7 ~ 1.9, respectively. 
As the genus increases, the invariant surface gets more 
and more complicated. Therefore, besides the integrable 
cases, the third special case which has a more "simpler" 
invariant surface, is expected for the pseudointegrable 
case 7 ~ 1.9 (we do not use 7 ~ 0.1 because we are 
interested in values of 7 in the interval [1,4]). Later we 
will come back to the special values 7 = 1.0, 1.9,3.0. 

The momentum distribution looks quite different for 
irrational mass ratios. Although the motion can be cr- 
godic, the number of momenta increases very slowly with 
longer system evolution 22] . However, this aspect is not 
in our present focus. 



III. TWO PARTICLES IN A ID-BOX WITH 
YUKAWA INTERACTION 

It is adequate to introduce the center-of-mass and rel- 
ative coordinates 



R = 



miqi + m2q2 
M ' 



and r = qi - q2, 



(2) 



respectively, with the total mass A/ — mi + m2 and the 
reduced mass ^ = mi 7712/ (mi -1-7712). In these new coordi- 
nates, the equations of motion describe a single compos- 
ite particle in the hyperspace (r, R) , called a hyperbilliard 
[rj. The Hamiltonian in relative coordinates is given by 



cos (9) = 



1 



1 - m2/mi 
l-|-m2/mi 1 -f- 7 



= A, 



(1) 



More specifically, writing 9 = ^tt, where m and n are in- 
tegers, at most An distinct velocity values occur. On the 
other hand, when 9 is an irrational multiple of tt, the ve- 
locities become uniformly dense [IB] in velocity space. As 
a consequence, it is at least possible for the two-particle 
hard point gas to be ergodic in velocity space if 9/^ is 
irrational. Although Casati and Ford did not show 
explicit results for irrational multiples of tt, they argued 
that their numerical results provide evidence support- 
ing ergodic behavior for irrational ^/tt by demonstrating 
that an increasing number of velocities is observed for a 
sequence of rational 9/^ approaching an irrational 9/^- 
value. Note however, that for irrational 9/t: infinite time 
may be required to observe all velocities. 

Although there are infinitely many mass ratios which 
give rational values of 9/tt, some of them are special. 
First, the integrable cases [13] 7 = 1,3 (or 1/3), which 
have 9 = and 9 = (or 7r/3), respectively. Relating 
Eq. ([T|) with results for the triangle billiard [I^] (or even 
to the polygonal billiard 0, [13, [lH ) it is possible to show 
that for the integrable cases the genus is equal g = 1 (the 
invariant surface of the billiard flow is a torus). For all 
other rational 9/Tr the dynamics is pseudointegrable [T3 ] 
and the invariant flow is not a torus (1 < 5 < 00). For 



H = 



2M 



El 
2^1 



(3) 



where the Yukawa potential V{r) has strength Vq and the 
parameter a > gives the interaction range vq = 1/a. 
Using scaled coordinates defined by (a 0) 

r = Tor, p = p^/E, R = roR, 



p = PVe, dt 



ro 



dr. 



and dividing Eq. ([3]) by E, the scaled new Hamiltonian 
is 



with 



p2 ~2 ^ 

H = + — + V(f) = e = 1, 

2M 2)1 ^ ' 



V{f)^V„^—, and Vb-^, 
r roE 



(4) 



(5) 



and scaled energy e = 1. When a = (tq —f 00) the 
transformation is independent of rg and Vq = ■ Under 
V{f) the composite's particle relative motion is subject 
to the force Q{f) = —dV/dr while its center of mass 
motion in R is free. For the case of equal masses the 
chaotic motion of ^ was already analyzed [1]. 
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A. Dynamics in tangent space, Lyapunov 
exponents (LE), and the origin of chaotic motion 

This section is dedicated to the analytical calculation 
of the Lyapunov spectrum. LEs are very useful! to de- 
scribe the dynamics in complex systems [l^l- When the 
motion is chaotic, at least one LE is positive. Its value 
is determined through the dynamics in tangent space, as 
will be shown below. 

Between collisions with the walls, equations of motion 
have the form 



F(7) = (f,/?,^,y)* = («,y,Q(r),0)*, 



(6) 



and it is easy to see that center of mass momentum AdR 
is a constant of motion. In relative coordinates, the com- 
posite particle moves under the influence of the force 
Q{f). This is a one-dimensional motion which is reg- 
ular and integrable. Collision with left and right walls 
cause a breaking of the translational symmetry of the 
system, and as a consequence, center of mass momentum 
is not a constant of motion anymore. The effect of left 
(right) wall collisions lead to the following change in the 
phase space point 7^ = (P, 7?*, w', V^') before the collision 
to 7/ = {r-f ,R-^ ,V^) after the coUision, 7/ = Tikji, 
with 



/I 














1 














-A 




Vo 





9 M 

M 


A 



(7) 



where A = (mi — m2)/M. The label k = 1(2) is used 
for particle 1(2). Since the two particles can never pass 
each other in the one-dimensional billiard, it is assumed 
without loss of generality that particle 1(2) never collides 
with the right (left) wall. The complete time evolution in 
phase space can be formulated by integrating the equa- 
tion of motion between collisions with the walls and by 
taking into account 7^ from Eq. ([7]) each time the com- 
posite particle collides with the walls. 

In order to calculate Lyapunov exponents, the time 
evolution of an infinitesimal path difference (the nearby 
trajectory) in the scaled tangent space {S^) must be de- 
termined, given by 



Sj{t) = M{t)Sj{to), 
with the scaled monodromy matrix 

dj{T) 



M(t) 



(8) 



(9) 



Lyapunov exponents arc the average rates of growth or 
shrinkage of such infinitesimal changes that are the eigen- 
vectors of M, 



A,; 



lim 

T — ^00 



logM^(T) 



(10) 



where fiiir) is the i-th. eigenvalue of M. The matrix M 
can be written itself as a product of matrices for small 
time steps. Since the motion between collisions is regular, 
Lyapunov exponents are zero. The situation is different 
for collisions of the composite particle with the walls. 
We follow the algorithm developed by DcUago et at [l^l 
to formulate the equations in the scaled tangent space 
according to 



0% 



rlC 

F{^,)-F{Cm) 



(11) 



where C = T>k^i gives the transformation at collisions 
with the walls with Dj. from Eq. Stj is the delay 
in the collision time of the (infinitesimal) nearby trajec- 
tory with respect to the collision time of the reference 
trajectory. Between the collision time of the main tra- 
jectory and the collision time of the nearby trajectory, 
the composite particle moves under the influence of the 
force Q{r) and the delay time can be determined from 



dSk 
dE 



nik 



Ik 



1% 



Ik 



^, (12) 
Qk 



where Sk is the action of the k particle, is the posi- 
tion of the nearby trajectory when the main trajectory 
collides with the wall, and is the collision point of 
the nearby trajectory. The energy dependence for qk, 

qi ~ y^2mi(£' — \m2q2 — is obtained from the 

energy conservation. Using Eqs. ([2]) we have 



9i(2) — R 



M 



M 



(13) 
(14) 



and Eq. can be written as 



MdR 



Afl^-(-)i(2)m2(i)« 
m2(i)dr 



1(2) 



Ml/ - (-)i(2)m2(i)z;' 



(15) 



Since V and v does not depend on R, the integral in 
R can be determined analytically. After the integration 
of Eq. terms proportional to SR^ = i?o — Rn and 

Sr'' ~ tq — rn will appear (quadratic terms {SR)^ and 
((5r')^ can be neglected in the linear analysis) and it can 
be written in scaled relative coordinates as 



- Stj =tk = -^= AkSR' + BkSr\ (16) 



E 



4 



where 



^ ' A/y»-- (-)i(2)m2(i)W*' 
_ ^ _(_)i(2) rf. m2(i)df 

'^'^ Sr' Jfo My-(-)i(2)m2(i)^' 

Finally, using Eqs. (El),© and (fTBl) . in pT|) . the collision 
of the composite particle in tangent space with left (fc ~ 
1) and right walls (fc = 2) under a force (5(^) ~ —dV/dr 
takes the form 6jf = MkSji with 



Mfc = (-) 



/ 


-A 


2 









9 P 

M 


A 












-AkQAk 


-A 


2 


V- 






9 


A,/ 



.(17) 

where Afc = —[A + (—1)'^]. The determinant of Mfc is 
equal to 1 and eigenvalues are also equal to |1|. However, 
the matrix elements proportional to Ak,Q,Ak and Bk 
generate positive LE when the global monodromy ma- 
trix M is constructed. Therefore, the presence of the 
interaction force Q and terms and Bk, related to the 
time delay in the tangent-space collision dynamics, are 
essential for the chaotic properties of the system. Other- 
wise, if Q = 0, no positive LE can be obtained. 

Two limiting situations can be discussed, a) The long 
range interaction (a = 0.0, tq —^ oo): at each collision 
with the wall, the interaction force Q(r) is finite and pos- 
itive LEs are generated; b) The short range interaction 
(ro ^1): in general Q{r) is close to zero if particles are 
sufficiently separated and the LEs are zero. However, the 
dynamics becomes chaotic due to a double collision pro- 
cess. For example, assuming that particles 1 and 2 are 
moving close together towards the left wall with a mutual 
repulsion close to zero due to the short range nature of 
the interaction. As particle 1 collides with the left wall, it 
changes its direction and moves towards particle 2, inter- 
acting with it. If such a double collision occurs infinitely 
close to the wall, the terms QB^ and QAk are not neces- 
sarily zero, a positive LE is generated and chaotic motion 
appears. This crucial role of double collisions to generate 
chaotic motion was also observed in another model of two 
interacting particles . The chaotic motion induced by 
the soft potential will be discussed in more detail in the 
next section. 



B. Signatures of regular motion in the distribution 
of largest finite-time Lyapunov exponents 

From the description of the dynamics by the mon- 
odromy matrices as constructed in the last section, it is 
clear that chaotic dynamics is generated in the presence 



of a soft interaction potential. The interesting question, 
however, is, if any signatures of non-ergodicity from the 
hard-point collision, or stickiness, can be identified. To 
this end we have investigated the distribution P(At, 7) of 
the finite-time largest Lyapunov exponents [2^ At as a 
function of the mass ratio 7. In general, for infinite time, 
the LEs Atx3 are well defined and do not depend on initial 
conditions. This holds also true for reasonably large finite 
times, if the motion is ergodic and the Lyapunov spec- 
trum has good convergence properties. In quasi- regular 
systems, however, where the chaotic trajectory may ap- 
proach a regular island and can be trapped there for a 
while, the value of the local LE can decrease. This will 
affect the convergence of At which depends now on the 
initial conditions. On the other hand, it implies that the 
distribution P(At), calculated over many initial condi- 
tions, contains information about the amount of regular 
motion (and sticky trajectories) in phase space. Usually, 
for fully chaotic systems P{At,j) has a Gaussian distri- 
bution (see for example [26j and references therein). 




FIG. 1: Mean value of the finite-time largest Lyapunov expo- 
nent calculated over 400 trajectories up to time t = 10* and at 
scaled energy e = 1, for (a) long interaction range (ro 00), 
i. e. Vo = 0.1, and (b) short range interaction (ro — 0.1), 
i. e. Vo = 1.0, with the hard walls located at g = ±1. For 
each trajectory the largest LE is evaluated over 10'' initial 
conditions. 

The mean (At), shown in Fig. [ll(a)-(b), decreases 
monotonically from roughly 0.9 to 0.54 in Fig. [TJa) 
(ro 00), and from 0.80 to 0.53 in Fig. [^b) (ro = 0.1), 
over a change of 7 from 1.0 to 4.0. This means that the 
dynamics is getting more and more regular as expected 
since 7 — > 00 constitutes an integrable limit with the 
heavy particle at rest. Figure [2] shows the finite-time dis- 
tribution of the largest LE. P(At, 7), for the case of long 
range interaction (ro — > 00). It reveals two indications 
for increasingly regular motion with growing mass ratio: 
a) the value of (At (7)) itself decreases and b) an increas- 
ing number of initial conditions lead to At (7) close to 
zero (related to sticky trajectories), or converge exactly 
to zero (regular islands). At 7 = 4.0 for example, about 
15% of the initial conditions lead to At = 0. The gray 
points below the main curve are related to chaotic tra- 
jectories which were trapped for a while close to regular 
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islands. Some examples will be given below. An interest- 
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FIG. 2: (color online). Finite-time distribution of the largest 
Lyapunov exponent P(At,7) calculated over 400 trajectories 
up to time t — lO'* and for ro oo. With increasing P{At, 7) 
the color changes from light to dark (white over yellow and 
blue to black). 

ing feature in Fig.[2]is the change of the width of P(At, 7) 
around the most probable defined through 



5P(A,,7) 



dAt 



= 0. 



(18) 



For mass ratio between 7 ^ 1.5 and 7 ^ 2.2, many initial 
condition lead to the same A^. In this region, Af (7) has 
a maximum as a function of 7. In other words, almost 
all initial conditions converge to the same LE and the 
dynamics should approach an "ergodic-like" motion. In 
fact, in this region the gray points (related to sticky or 
trapped trajectories) below the main curve almost disap- 
pear. 

A similar behavior is found in the short interaction 
limit (ro = 0.1) shown in Fig. [21 Compared to Fig. [21 two 
main differences can be observed: firstly, the maximum 
for Af (7) from Fig. [21 (between 7 ^ 1.5 and 7 ^ 2.2) 
is divided into two maxima, one close to 1.5, and the 
other one close to 2.4. Therefore, a minimum of A((7) 
appears in between (7 ~ 1.9), where trapped trajectories 
are expected to occur more often if compared with the 
two maxima 7 ^ 1.5 and 2.4. Secondly, the abrupt ap- 
pearance of many gray points below the main curve close 
to 7 ~ 2.7; which may indicate that a regular island is 
born. This will be shown in more details below. 

A more systematic way to uncover this trend is to fol- 
low P(Af,7) = Pa (7) as a function of the mass ratio 7 
shown in Fig. [H (top). If Pa is large, a large fraction of 
initial conditions lead to the same At and trapped tra- 
jectories are rare. For example, the maximum of Pa (7) 
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FIG. 3: (color online). Finite-time distribution of the largest 
Lyapunov exponent P(At,7) calculated over 400 trajectories 
up to time t = lO** and for ro = 0.1. With increasing P(At, 7) 
the color changes from light to dark (white over yellow and 
blue to black). 



in Fig. [4] (top, ro — > 00, black curve) close to 7 ^ 1.9, 
is the region in Fig. [21 where gray points below the main 
curve are rare. The fast variation of Pa (7) is due to sta- 
tistical fluctuations in the determination of At over the 
400 initial conditions. However, two main valleys can be 
identified in the black curve of Fig. [H one close to 7 = 1 
and the other close to 7 = 3. These are exactly the mass 
ratios for which the hard point-like collision dynamics 
(scction[lI| is integrable. The gray curve of Fig.[4]for the 
short interaction (rp = 0.1), which is closer to the limit 
of hard-point collisions, presents an additional little val- 
ley at 7 ~ 1.9. This is the pseudointegrable case with 
genus g = 2 which appears in the hard-point collision. In 
other words, if 7 is close to values for which the dynam- 
ics in the point- like gas is integrable {g = 1) or "simpler" 
{g = 2), then the dispersion around Af increases so that 
Pa (7) decreases, and signatures of non-ergodicity are ex- 
pected under additional Yukawa interaction. Another 
interesting observation is that the minimum at 7 ~ 1.9 
(see gray curve from Fig. [4]) disappears in the long in- 
teraction limit ro ^ 00 (see black curve from Fig. [J]). It 
means that the regular motion from the integrable cases 
of the hard-point collision survives longer under the per- 
turbation of the soft interaction than the regular motion 
from the pseudointegrable case. 

For fully chaotic systems the quantity Pa(7) is just 
the maximum of a Gaussian distribution and it should 
increases linearly with t, since the variance a = ((Aj — 
(At)^))^/^ for such systems goes with 1/t. This behavior 
of cr has been observed by studying ergodicity in high- 
dimensional sympletic maps [27[, and used for the de- 
tection of small islands in the standard map (2^. For 
our case we found that the time dependence of Pa (7) 
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FIG. 4: Top: Normalized distribution Pa (7) of the most prob- 
able Lyapunov exponent Aj for ro — > 00 (black) and ro = 0.1 
(gray). Bottom: a for a) ro — > 00 and b) ro = 0.1 



goes with ^ t*^'^, which means that a significant number 
of islands is present. In order to compare a with Pa(7), 
Figs. [3^) and b) show the behavior of cr as a function of 7. 
Each time trapped trajectories are present, a should has 
a maximum, exactly the opposite behavior from Pa (7). 
In Fig.[4]a) we see two maxima, one close to 7 = 1.0 which 
corresponds to the minima of Pa (7) for the same value of 
7, and the other one at 7 3.5 which as no analogous for 
Pa (7). For Fig.|3]b) we again sec two maxima, one close 
to 7 = 1.0 and the other one located between 7 ~ 3.0 
and 7 ~ 3.5. Besides the peaks at 7 = 1.0, the other 
peaks are not located at the integrable (7 = 3.0) and 
pseudo-integrable (7 ~ 1.9) cases from the hard point- 
like collisions cases. The peak around 7 ~ 3.0 in Fig. 2] 
b) is correct, but by far not precise. This shows that the 
quantity a is not sensible enough to detect signatures of 
the regular dynamics from the hard-point collision case. 
The reason for this is that a is not able to detect small 
regular islands which appear in the phase space, as for ex- 
ample the points related to zero LEs in Fig.[3]for 7 = 1.9. 



It is worth to mention that the maxima and minima from 
^a(7) does not change significantly with the number of 
initial conditions. 

We can conclude that signatures from regular struc- 
tures in phase-space exist and are uncovered by the dis- 
persion of the largest Lyapunov exponent At most clearly 
visible along the cut Pa (7) defined by the number of oc- 
currences of the most probable Lyapunov exponent 
(Eg. [T8)) . No such signatures, however, are visible in less 
sensitive quantities such as the mean Lyapunov exponent 
(Af) or the fluctuation cr. 




FIG. 5: Poincare Surfaces of Section for ro — > 00 and a) 7 = 
3.0 and c) 7 = 4.0. Figures c) and d) are the corresponding 
magnifications of the quasi-regular regions. 

The interpretation of the Lyapunov properties are sup- 
ported by the relevant PSS. Figured] (top) shows PSSs for 
interaction rg — > 00 and cases 7 = 3.0 and 7 = 4.0, where 
Pa (7) has a minimum and maximum (see black curve in 
Fig.[4|). respectively Although the system is more chaotic 
for 7 = 3.0 than for 7 ~ 4.0, trapped trajectories appear 
near the island for 7 = 3.0, see the corresponding mag- 
nification in Fig. [5] (top, right). Such trapped motion 
near regular islands, which does not appear for 7 = 4.0 
(see Fig. [Sl^bottom, right)), affects P(Af,7) and conse- 
quently. Pa (7) has a minimum near 7 ~ 3.0. Another 
example is shown in Fig. [S] for tq — 0.1. For 7 = 1.5 
(top, left), which has a maximum in Fig. 3] (top, gray 
curve), no trapped trajectories are found up to the time 
propagated. For 7 = 1.8 (top, right) which has a min- 
imum in Fig. |4] (top, gray curve), trapped trajectories 
start to appear. Moreover, the abrupt appearance of gray 
points for 7 ~ 2.7 below the main curve Fig. [3] can be 
nicely explained using the PSS. Figure compares the 
PSSs for 7 = 2.6 (bottom, left) with 7 = 2.8 (bottom, 
right), where gray points appear in Fig. [31 Clearly, it 
can be seen that when a regular island is born trapped 
trajectories start to appear around the island and, as a 
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FIG. 6: Magnification of tiie Poincare Surfaces of Section with 
ro = 0.1 for a) 7 = 1.5, b) 7 = 1.8, c) 7 = 2.6 and d) 7 = 2.8. 



consequence, many initial conditions with lower LE are 
obtained (gray points in Fig. [4]). The regular island is 
related to a period-4 orbit which has the following prop- 
erty: for each second hit of particle 1 with the left wall, 
particle 2 is at rest. There is a similar periodic orbit for 
the hard-point collision case [20| . 

Note that although Pa (7) is determined only from tra- 
jectories with positive At, it provides information about 
the amount of regular structure in phase space through 
chaotic trajectories which are trapped close to regular 
islands. Since trapped trajectories are characteristic for 
mixing in phase space, Pa (7) provides a tool to analyze 
phase space mixing. 



IV. TWO PARTICLES IN A CIRCULAR 
BILLIARD WITH YUKAWA INTERACTION 

In order to show the utility of Pa (7) for systems with 
higher dimensions, we discuss now the case of two in- 
teracting particles in a circular billiard. The interaction 
between particles is still of Yukawa type. The chaotic 
motion is now generated by the combined effect of the 
curvature of the walls from the circular billiard and the 
double collisions discussed in last section. The phase 
space is 8-dimensional and it is not possible to construct 
an adequate PSS which allows to look at the underline 
dynamics. Trajectories will fill out any chosen PSS and 
no information about details of the dynamics can be ob- 
tained. Sticky trajectories, for example, which may cause 
non-ergodicity due to a partial focusing of trajectories 
[sot . 3.re difficult to detect. In such partial focusing, an 



infinitesimal family of nearby trajectories that starts out 
parallel will lead to Lyapunov exponents which converge 
very slowly in time. This is a typical behavior in high di- 
mensional quasi-regular systems. In this section we show 
the effectiveness of Pa (7) to obtain relevant informations 
in high-dimensional quasi-regular systems. 

Figure [7^) shows the finite-time distribution of the 
largest LE for the case of long range interaction (rg — > 00) 
in the circular billiard. The value of the the mean LE cal- 
culated of the 400 initial conditions, decreases growing 
the mass ratio and the regular motion increases. This 
is the only information we can get from the LE about 
the complicated dynamics of the two interacting parti- 
cles inside the billiard. However, the gray points below 
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FIG. 7: (color online), a) Finite-time distribution of the 
largest Lyapunov exponent P(At,7) calculated over 400 tra- 
jectories up to time t = 10* and for ro — > 00 for the circular 
billiard. With increasing P(At, 7) the color changes from light 
to dark (white over yellow and blue to black); b) Normalized 
distribution Pa (7) of the number of occurrencies of the most 
probable Lyapunov exponent . 



the main curve of Fig. [TK) are related to chaotic tra- 
jectories which were trapped for a while close to regu- 
lar islands. Therefore some additional information about 
sticky trajectories and the amount of quasi-regular mo- 
tion in phase space may be obtained from Pa (7). This 
is shown in Fig. EId). One minima is observed close to 
7 ~ 1.0 where trapped trajectories and partial focusing 
of trajectories in phase space are expected. On the other 
hand, a more ergodic-like motion is expected for 7 ~ 3.0, 
where Pa (7) has a maximum. This is different from the 
results for the ID-Box from last section (see the mini- 
mum for r — *■ 00 at 7 = 3.0 in Fig. [3]Top). In other 
words, we clearly see that signatures from the integrable 
case 7 = 3.0 from the ID hard point-like collision case 
disappear in the two dimensional case. This tells us that 
the dynamics from the integrable case 7 = 3.0 is a con- 
sequence that particles cannot pass each other in the ID 
case. As the masses ratio increase to 7 = 4.0, the LE 
decreases and the system is more regular. In this case 
Pa (7) also decreases slowly and the amount of regular 
islands, and the motion around them, increases. 
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V. CONCLUSIONS 

Usually, chaotic motion is generated through non- 
linear equations of motion. As a consequence, exponen- 
tial divergence of nearby trajectories occurs which can 
be quantified by positive Lyapunov exponents. Another 
source of chaotic motion emerges even for linear equa- 
tions of motion through boundary conditions. The ad- 
vantage of such systems is the possibility to obtain rigor- 
ous mathematical results. One example is a point parti- 
cle moving among high-dimensional cylindrical scatterers 
[2^ . which is similar to the high-dimensional Lorentz gas. 
In these systems, strongly chaotic motion is generated 
due to the convex curvature of hard disks or spheres. In 
fact, the collision time delay between nearby trajectories 
due to the curvature of the surface, is responsible for the 
chaotic motion. 

Curved boundaries are not present in the one- 
dimensional confinement considered in this paper. For 
the hard point-like collision case non-ergodic motion is 
generated when the mass ratio gives a 9 value [from 
Eq. ([1])] which is a rational multiple of tt. Ergodic mo- 
tion, on the other hand, may be obtained for irrational 
multiples of O/tt. Wc have shown that for additional soft 
Yukawa interaction between the two particles, chaotic 
motion is obtained for any mass ratio. Double collisions 
of particles, which occur very close to the walls, are es- 
sential to generate positive LEs in the short interaction 
range limit. The coUisional time delay in tangent space, 
together with the soft Yukawa interaction, are responsi- 
ble for the chaotic motion. The mean of the largest finite- 
time LE (At) decreases smoothly as the mass ratio in- 
creases and does not provide detailed information about 
the phase space structure. This type of information is 
provided by the probability distribution of the largest 
finite-time Lyapunov exponent P(At,7). It reveals that 
the dispersion around (Aj) increases when trapped tra- 
jectories arc present in the phase space. Wc have shown 
that a cut through P(Ai,7) along the number of occur- 
rencies of the most probable Lyapunov exponent. Pa (7), 
gives a quantitative measure of the influence of regular 
motion in mixed phase space. Specifically for the sys- 
tem studied here, we have shown that Pa (7) decreases 



when the structure in phase-space is more regular and the 
mass ratio is close to the integrable cases 7 = 1,3 (genus 
g = 1) or to the "simpler" dynamics (pseudointegrable 
at 7 1.9, with g = 2) from the hard point-like collision. 
We also observed that the regular motion of the inte- 
grable cases from the hard-point collision survives longer 
under the perturbation of the soft interaction than the 
regular motion from the pseudointegrable case. Hence, 
the dynamics under the additional Yukawa interaction, 
although in principle chaotic, "remembers" the integrable 
and pseudointegrable dynamics in the system without the 
soft Yukawa interaction. This is certainly a subtle effect 
and therefore, we expect that in general the number of 
occurrencies of the most probable Lyapunov exponent 
provides a sensitive tool to probe details in phase space 
dynamics. We have also shown that this quantity is much 
more sensitive if compared with the mean square fluctua- 
tions of the LE. In contrast to Poincare sections this tool 
is easily applicable in higher dimensional systems, where 
trapped trajectories may cause non-ergodicity due to a 
partial focusing of trajectories [10]. In order to show this 
we calculated P(A(,7) for two interacting particles in a 
circular billiard, where the phase space is 8-dimensional 
and it is not possible the construct an adequate Poincare 
section to analyze the dynamics. We show that by in- 
creasing the mass ratio, the mean LE decreases and the 
system gets more regular. Furthermore, for a minimum of 
Pa (7) at 7 ^ 1.0 trapped trajectories are expected, and 
a maximum at 7 '--^ 3.0 ergodic-like motion is expected. 
Therefore, Pa (7) can be used in higher-dimensional sys- 
tems as a tool to characterize the dynamics. 
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